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文章 编号 :1005-3085(2010)05-0827-06 
求 线性 矩阵 方程 双 对 称 最 小 二 乘 解 的 变形 共 斩 梯 度 法 


田 小 红 ， 张 凯 院 
(西北 工业 大 学 应 用 数学 系 ， 西 安 710072) 

摘 要 : 本 文 基于 求 线性 代数 方程 组 的 共 辆 梯度 法 的 思想 ， 通 过 特殊 的 变形 与 近似 处 理 ， 建 立 了 求 一 般 线 
性 矩阵 方程 的 双 对 称 最 小 二 乘 解 的 迭代 算法 ， 并 证 明了 迭代 算法 的 收敛 性 。 不 考虑 舍 入 误差 时 ， 
迭代 算法 能 够 在 有 限 步 计算 之 后 得 到 矩阵 方程 的 双 对 称 最 小 二 乘 解 ， 选取 特殊 的 初始 矩阵 时 ， 还 
能 够 求 得 矩阵 方程 的 极 小 范 数 双 对 称 最 小 二 乘 解 。 同 时 ， 也 能 够 给 出 指定 矩阵 的 最 佳 逼 近 双 对 称 
和 矩阵。 算 例 表明 ， 和 迭代 算法 是 有 效 的 。 
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1 引言 


和 矩阵 方程 特殊 解 的 计算 问题 在 电学 、 结 构 动力 学 、 振 动 理 论 、 自 动 控制 理论 等 领域 都 有 重 
要 应 用 。 文 献 [1-6] 建立 了 求 矩 阵 方 程 AXB = CRE AXB +CXD = 下 某 些 特殊 解 的 迭代 方 
法 。 本 文 建 立 求 一 般 线 性 矩阵 方程 双 对 称 最 小 二 乘 解 的 变形 共 斩 梯 度 法 (MCG 算法 )。 

R” RI mxn KERERE, SR” Aog nr SOWMEOBEESE S. AG B RIPER A 5 BI 
Kronecker 积 ， 定 义 实 和 矩阵 4 与 BB 的 内 积 为 [4, B] = trace(AT B), 由 此 导出 和 矩阵 的 Frobenius 范 
HJA = Vv[4, 4]，vec(4) 表 示 将 矩阵 4 按 行 拉 直 构成 的 列 向 量 ，e; 表示 单 位 矩阵 ,的 
第 i 列 。 记 S= (en,en_1,… ,e) WA S = n, ST = So 

定义 1 FEX c R” 的 元 素 满 足 


Tj = X4 = Tntl jntliy 1, J = 12,-,n, 


则 称 X ANIERE. E n MAAIKE PERI E rU BIR., 
XT —BRRUS £c TEXE EE 77 FE 


AıXBı +---+ ANXBN =C, 
研究 下 面 两 个 问题 : 
问题 1 设 A;€R™*", B,€ Rn"xP ((212,-.,N), Ce Rm*?, RX € BSR"*" 使 得 
N 
| > AXB -C| — min. (1) 
i=l 


问题 2 MEXER, Sp 表示 问题 1 的 解 集合 ， 求 总 € Ss 使 得 
IX -Xl = min |X- XI. (2) 
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2 问题 1 的 等 价 转化 
引 理 1 REX c BSR"*" 的 充 要 条 件 是 和 X= XT —SXS. 
引 理 2 ” 若 矩 阵 X€ SR”, WX--SXSc BSR^*", 
为 了 讨论 方便 引进 记号 
N N 
f(X)= 32M (474;XB;B] + BiB} XA} A; 
i=1 j=1 
+SATA;SXSB;B!S + SBiBISXSAT 4;5) 
N 
Q= M (A7 CB] + BiCTAi+SATCBIS + SBiC™ A,S) 
i=1 
定理 1 求解 问题 1 等 价 于 求 矩阵 方程 
f(X)- (3) 
的 双 对 称 解 ， 而 且 该 矩阵 方程 一 定 有 双 对 称 解 


证 明 C4 X e BSRnxn 时 ， 由 引 理 1 知 X = XT = SXS。 因 此 ， 求 X e BSR"x" 使 
得 (1 成立， 等 价 于 求 Xe BSR"*" 使 得 
N 2 N 2 
| xxn -cl +Z arxar - 7] 
gs] i—1 


N 2 N ; 
+|D 4sxsB: - c| «| Y 87sx847 - c7 | = min 
5 i=1 


(4) 
下 面 证 明 求 极 小 值 问题 (4) 的 双 对 称 解 等 价 于 求 矩阵 方程 (3) 的 双 对 称 解 。 将 矩阵 方程 组 


Mz 


A:iX B; =C, 


e. 


m 


x BTXAT = CT, 


e. 


(5) 
A;SXSB, = C, $ BISXSAT = CT, 
1—1 i—l 
按 行 拉 直 可 得 线性 方程 组 
N 
X 4:8 BT 
1—1 
N vec(C) 
IO nA sec(C7) 
E - vec( X) d ; (6) 
E (4:8) e (B7 S) EAG) 
Ps ve«CT) 
XS) 8 (4:5) 
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易 见 ， 求 线性 方程 组 (6) 的 最 小 二 乘 解 ， 等 价 于 求 矩 阵 方程 组 (5) 的 最 小 二 乘 解 ， 也 就 是 求 极 小 
值 问题 (4) 的 解 。 记 线性 方程 组 (6) 的 系数 矩阵 为 了 ， 那 么 


N N N " 
vT = (3 AT & B, 3 Bio AT, Y SAT) 8 (SB), Y (5B) 9 (SAD)). 
fet Scl i=1 i=l 
引进 记号 
N N 
w = Y Y ((474;) € (BBT) + (B:BT) e (AT Aj) 


+ (SAT A;S) 8 (SB;:B7 S) + (SB:B7 8) & (SAT A;8)), 
线性 方程 组 (6) 的 正规 方程 组 可 写 为 
W -vec(X) = V" ( (vec(C))* , (aec(C7))T， (vez(C))7, (vee(CT™))T), (7) 


将 正规 方程 组 (7) 还 原 为 矩阵 方程 就 是 矩阵 方程 (3)。 
下 面 证 明和 矩阵 方程 (3) 有 双 对 称 解 。 因 为 矩阵 方程 (3) 有 解 ， 所 以 可 设立 是 它 的 一 个 解 (未 
必 是 双 对 称 解 )， 那 么 
f(X) - Q, (8) 
4 
g(X) = Tx (Ñ+ X7 + St + RT)S), 


由 引 理 2 知 g(X) e BRS”, RB) 取 转 置 ， 然 后 分 别 左 、 右 乘 S， 联 立 运算 可 得 f(g(X)) = 
Q， 这 表明 o (X) 是 矩阵 方程 (3) 的 一 个 双 对 称 解 。 


3 ”求解 问题 1 8935 TE CS Bb HE 


下 面 建立 求 矩 阵 方程 (3) 双 对 称 解 的 MCG 算法 ， 也 就 是 求解 问题 1 的 MCG 算法 。 

步骤 1 "BEA; € R”, B; eR?” (i=1,2,--,N), Ce Rm™*?, Xo e BSR"*^, 
步骤 2 WH Ro =Q- f(Xo), Go = f(Ro), Po = Go. fifi k := 0; 

步骤 3 WRR =0, Bib FW, Vk:—kocid; 

步骤 4 计算 

l4 i? 


Xy = Xk-1 + 
Eeki Ppa? 





Pk-1, 


| il? 


TRUSTe a 





f (Pe-1), 





下 Pe, (变形 之 一 ) 


步骤 5 返回 步骤 3。 

引 理 3 X,Y e R?x*?， 则 有 [fF(X),Y] = [X, fY). 

定理 2 ”如果 X* 是 矩阵 方程 (3) 的 双 对 称 解 ， 那 么 对 任意 初始 双 对 称 和 矩阵 X9, MCG 算法 
PRIER Xa, Ri 和 Gi 满足 [Gi, X* — Xi] = |R: (i= 0,12….)。 
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证 明 因为 X* 是 矩阵 方程 (3) 的 双 对 称 解 ， 所 以 f(X*) = Q@， 当 i = 0,1,2,:… 时 ， 有 
[Gi, X* - Xi] = [f(), X* — Xi] = [R;, f(X* — X] 
= [R5 f(X*) - f(X)] = [Ri, Q ~ f(X))] = [Ri, Ri] = l4. 


定理 3 如 果 X* 是 矩阵 方程 (3) 的 双 对 称 解 ， 那 么 对 任意 初始 双 对 称 和 矩阵 Xo，MCG 算法 
"X; Ri M P, WE [P X* — Xj] = [RJ (i=0,1,2,.…)。 

利用 定理 2 采用 数学 归纳 法 即 可 证 得 结论 。 

推论 1 X MCG FAF RERE R RIP; 4$ Ri 40, WP; #0(i= 0,1,2, )o 

定理 4 对 MCG 算 法 中 的 矩阵 Ri 和 PP， 车 有 正 数 k， 使 得 R; Z 0(i = 0,1,:… ,k)， 
则 [Ri, R;] =0, [Pi, P] =0(i Æj; i,j 90, 1,--- , k). 

证 明 根据 矩阵 内 积 运算 的 交换 律 ， 只 需 对 0 < i < ; < 证 明 结 论 成 立即 可 。 

SERI 采用 数学 归纳 法 证 明 [Ri, Riy] = 0, (P, PH = 0(i= 0,1,- k)s 

SEA 假设 当 0 <i< 和 有 111<1I « kB, [Ri Rinu] = 0，[ 忆 ,PH = 0 成 立 ， 证 
88 [Ri, Riy] = 0, [Pi, Piti] = 0。 详 细 过 程 略 。 

定理 5 对 任意 初始 双 对 称 矩 阵 Xo， 和 矩阵 方程 (3) 的 双 对 称 解 可 在 有 限 步 计算 后 得 到 。 

这 是 因为 在 有 限 维 矩阵 空间 Rmx 中 ，Ro, Ri, Ro, 两 两 正 交 ， 所 以 必定 存在 正 整数 k < 
2， 使 得 Ri = 0。 由 于 实际 计算 过 程 中 不 可 避免 的 舍 入 误差 影响 ，MCG 算法 并 不 能 在 有 限 步 
迭代 计算 后 得 到 问题 1 的 精确 解 ， 所 以 只 能 把 它 作 为 迭代 算法 使 用 。 

3147 RA e R", be Rm， 线 性 方程 组 hz = 2 有 解 ， 则 极 小 范 数 解 z* 唯一 ， 
且 z* = Atb € R(AT), WERT) PRA Az = 8 的 一 个 解 。 这 里 4+ 表示 4 的 Moore- 
Penrose 3f 


定理 6 EEDI (3) AEAII, F ARDINI ERE 
Xo = f(H), He BSR"™Y", (9) 


则 MCG FA EA RhE BEREA TE (3) 的 极 小 范 数 双 对 称 解 。 
证 明 ”将 矩阵 方程 (3) 按 行 拉 直 可 得 线性 方程 组 


W .vec(X) = vec(Q), (10) 


按照 式 (9) 选取 初始 双 对 称 和 矩阵 Xo 时 ， 由 MCG 算法 得 到 的 矩阵 Xn 具有 类 似 于 Xo 的 表示 形 


Ved( Xk) = W -vec(Y,) e RW) = R(W™), 
由 定理 5 知 ， 存 在 加 使 得 矩阵 方程 (3) 的 双 对 称 解 X* = Xr。， 从 而 vec(X*) 是 线性 方程 


组 (10) 的 解 。 由 Tec(X*) = vee(X4,) € R(WT) 及 引 理 4 可 得 vec(X*) = Wt .vec(Q@) 是 线 
性 方程 组 (10) 的 极 小 范 数 解 ， 从 而 X* 是 矩阵 方程 (3) 的 极 小 范 数 双 对 称 解 。 





4 问题 2 的 等 价 转化 


对 任意 给 定 的 矩阵 吉 e R"x"， 问 题 1 的 解 集合 Szg 非 空 ， 若 X e Se C BIR”, RE 
对 称 和 矩阵 与 反对 称 和 矩阵 的 正 交 性 ， 以 及 双 对 称 和 矩阵 (Y = YT = SY S) 与 对 称 次 反对 称 矩 
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BE (Y = YT = 一 SY5) 的 正 交 性 ， 可 得 








Ix - XIP 
pee IST 
| 
-ep ET 


于 是 ， 问 题 2 中 求 和 es Sg 使 


IŽ- ži pin IX- Xll 
ESE 














TMT R X € SE 使 
n Y YT Y yT y v1 X XT S 
|2 X+XT+S(X+X S| = min |x X +XT+S(X t X!) l, (11) 
4 X€Sg 4 
令 
xa 元 TT Ya xT N vao yT va XT 
X-x 天 十 -S(X--X E CC voa EE +S(X +X )S 5.. 
4 i=1 4 
4 X € Sg Ff, WEDE (3) 可 写 为 
N 
f(X)= > (ATCBT + BiC" A; + SATÓBI S + SBiC" A,S), (12) 
i=1 


矩阵 方程 (12) 与 矩阵 方程 (3) 的 差别 只 是 左 端 矩阵 不 同 。 根 据 式 (11). RA A2 E FRIE 

阵 方程 (12) 的 极 小 范 数 双 对 称 解 羡 。 取 初始 矩阵 头 o。 = f(H), HP HÀ BSR?” FERE 

阵 ， 使 用 MCG 算法 可 求 得 矩阵 方程 (12) 的 唯一 极 小 范 数 双 对 称 解 六 "， 从 而 得 到 问题 2 的 解 
ue E pe 


4 
算 例 结果 表明 该 MCG 算法 是 有 效 的 。 
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The Modified Conjugate Gradient Method for the Least Squares 
Bisymmetric Solution of the General Linear Matrix Equation 


TIAN Xiao-hong, ZHANG Kai-yuan 


(Department of Applied Mathematics, Northwestern Polytechnical University, Xi'an 710072) 


Abstract: In this paper, an iterative method is proposed to find the least squares bisymmetric solutions 
of the general linear matrix equation. This method is obtained through a special transformation and the 
approximate disposal on the base of conjugate gradient method for solving linear algebraic equations. 
The convergence of this method is also given. By using this method, a least squares bisymmetric 
solution can be obtained within the finite iterative steps in the absence of round-off errors, and the 
solution with least norm can be obtained by choosing a special initial bisymmetric matrix. In addition, 
its optimal approximation solution to a given matrix can be obtained. Numerical examples show that 
the iterative method is efficient. 

Keywords: bisymmetric matrix; least squares solution; least-norm solution; iterative method; optimal 


approximation 
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